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A CONJECTURE ON EXCEPTIONAL ORTHOGONAL POLYNOMIALS 



DAVID GOMEZ-ULLATE, NIKY KAMRAN, AND ROBERT MILSON 

Abstract. Exceptional orthogonal polynomial systems (X-OPS) arise as eigenfunctions of Sturm-Liouville 
problems and generalize in this sense the classical families of Hermite, Laguerre and Jacobi. They also 
generalize the family of GPRS orthogonal polynomials introduced by Gariiiena et al., |3]- We formulate 
the following conjecture: every exceptional orthogonal polynomial system is related to a classical system 
by a Darboux-Crum transformation. We give a proof of this conjecture for codimension 2 exceptional 
orthogonal polynomials (X2-OPS). As a by-product of this analysis, we prove a Bochner-type theo- 
rem classifying all possible X2-OPS. The classification includes all cases known to date plus some new 
examples of X2-Laguerre and X2-Jacobi polynomials. 



Q ' 1. Introduction 

cn ' 

The past several years have witnessed a considerable level of research activity in the area of exceptional 
orthogonal polynomials, which are new complete orthogonal polynomial systems arising as eigenfunctions 
of Sturm-Liouville operators, extending the classical families of Hermite, Laguerre and Jacobi. The first 
examples of exceptional orthogonal polynomial systems were discovered in [12' and [131 as a result of the 
T^ ', development of a direct approach [8^ to exact or quasi-exact solvability for spectral problems in quantum 

j^ ' mechanics that would go beyond the classical Lie algebraic formulations [T7 j [22 l l34] . 

The exceptional orthogonal polynomial systems and the Sturm-Liouville problems that define them 

have some key poperties that distinguish them from the classical orthogonal polynomial systems, and 

which we would like to highlight. The most apparent one is that they admit gaps in their degrees, in the 

sense that not all degrees are present in the sequence of polynomials that form a complete orthonormal 

set of the underlying weighted L^ space, even though they are defined by a Sturm-Liouville problem. This 

iy-s ' means in particular that they are not covered by the hypotheses of Bochner's celebrated theorem on the 

00 ! characterization of orthogonal polynomial systems defined by Sturm-Liouville problems [2] . 

^^0 ' The number of gaps in the sequence of degrees of the polynomials appearing in a complete family 

fT^ I will be referred to as the codimension and we will use the symbol X^ to denote the various complete 

^D ' orthogonal systems of codimension m. A second order differential operator is exceptional if it preserves 

^N I some exceptional polynomial fiag, but does not preserve the standard polynomial fiag generated by the 

monomials. Thus, and in contrast with the classical families where the defining differential operator 

has only polynomial coefficients, the second order differential operators corresponding to the exceptional 

k>( ' families have poles in their coefficients, although all their singular points happen to be regular. 

'j_j I The first explicit examples of families of exceptional orthogonal polynomials are the Xi-Jacobi and 

C^ ' Xi-Laguerre polynomials, which are of codimension one, and were first introduced in [12] and [l^. In 

these papers, a characterization theorem was proved for these orthogonal polynomial families, realizing 

them as the unique complete codimension one families defined by a Sturm-Liouville problem. One of the 

key steps in the proof was the determination of normal forms for the fiags of univariate polynomials of 

codimension one in the space of all such polynomials, and the determination of the second-order linear 

differential operators which preserve these flags pT | [T6 ] . 

It is Quesne |25l[^ who first observed the presence of a relationship between exceptional orthogonal 
polynomials, the Darboux transformation and shape invariant potentials. This enabled her to obtain 
examples of potentials corresponding to orthogonal polynomial families of codimension two, as well as 
explicit families of X2 polynomials. Higher-codimensional families were first obtained by Odake and 
Sasaki [28]. The same authors further showed the existence of two families of X^-Laguerre and X^- 
Jacobi polynomials, the existence of which was explained in [T3] for Xm-Laguerre polynomials and in [IB] 
for Xm- Jacobi polynomials, through the application of the isospectral algebraic Darboux transformation 
first introduced in [9l[T0]. We also refer to [33] for similar results, and to [14l[T6] for the proof of the 
completeness of the X^-Laguerre families. We also note that some examples of exceptional Hermite 
polynomials were known to the quantum physics community in the early 90s, [4j, and are actively studied 
today under the name of CPRS systems, [3j[7]. It should as well be noted that the exceptional Laguerre 
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polynomials have already been used in a number of interesting physical contexts, for Dirac operators 
minimally coupled to external fields, |21j . or in quantum information theory, [BJ. 

The papers cited above contain many examples of orthogonal polynomial families of arbitrary codi- 
mension arising from the Laguerre and Jacobi system by the application of the Darboux transformation. 
However, as was shown in [TF, the above list is not exhaustive: novel exceptional polynomials can be 
constructed by means of multi-step Darboux or Darboux-Crum transformations. The multi-step idea 
was further applied to exactly solvable and shape-invariant potentials up in |181I27|[52] . However, an 
essential question that remains open is to know whether these families exhaust all the possibilities of 
higher-codimensional complete orthogonal polynomial systems, in other words whether all the higher- 
codimensional complete orthogonal polynomial systems are generated by the application of successive 
algebraic Darboux transformations. We conjecture this result to be true. In order to prove such a result, 
one approach would be to try to carry out for all codimensions an analysis similar to the one performed 
in |llfll3) in codimension one, identify the complete orthogonal sets amongst the resulting families and 
show that all of these can be obtained from the classical codimension-zero families by iterating algebraic 
Darboux transformations (we will refer to these as multi-step Darboux transformations). This seems like 
a very challenging task in the absence of a general classification strategy that would lead to normal forms 
for flags of univariate polynomials for all codimensions. Even in the codimension two case, the question 
would be quite difficult to answer if we were only using the tools that were at our disposal in [llj . We are 
nevertheless able to give a complete answer to this question for codimension two families by suitably refin- 
ing the approach taken in these earlier papers. In particular the possible pole structure of the coefhcients 
of the operators that preserve the codimension two fiags plays a key role in the analysis. 

Since the main objects of our study are orthogonal polynomial systems that arise as eigenfunctions of 
a Sturm-Liouville problem, let us give a definition: 

Definition 1.1. We define a Sturm-Liouville orthogonal polynomial system (SL-OPS) as a sequence of 
real polynomials yi (x) , ?/2 (x) , j/3 (x) , . . . , with deg yi > deg j/j if i > j , satisfying the following conditions: 

(i) There exists a second order differential operator 

T[y] = p{z)y" + q{z)y' + r{z)y 

with rational coefficients p,q,r such that T[yi\ = Xiyi for all i, with the Xi distinct, 
(ii) There exists an interval / ~ (a, b), — cx) < a < fe < 00 such that the weight function 

W(x) = —^^ exp I / -dx 
P{x) \J p , 

is positive, that is W{x) > for x (z I, such that all moments are finite, 

x'W{x)dx<oo, i = 0,1,2,3..., 

and such that p{x)W{x) — ?> at the endpoints x ^ a,b. 
(iii) The polynomial sequence is complete, meaning that span{yi}°Zi is dense in L'^{Wdx,I); 

Remark 1.1. It follows from the above definition that the operator T is essentially self-adjoint on the 
weighted Hilbert space L^{I,Wdx) and that the eigenpolynomials are orthogonal, meaning that 

l-b 



/ Wytyjdx = hSij, h > 0, 

J a 



for some constants fc,-. 



Remark 1.2. If degy^ = i — 1 for all i, we are dealing with one of the classical orthogonal polynomial 
systems of Hermite, Laguerre and Jacobi: the polynomials in question span the standard polynomial flag 
and p, g, r are polynomials of degrees 2, 1 and respectively, [5]. 

Remark 1.3. If the degree sequence {degyi}^]^ does not contain all non- negative integers, then we will 
have an exceptional polynomial system (X-OPS), and the coefficients of T will be purely rational (as 
opposed to polynomial) functions. 

Remark 1.4. We shall see in Section 5.1 that the eigenvalue equation T[y] = Xy can be put into Sturm- 
Liouville form. 
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Even though several famihes of X-OPS have now been described in the hterature, the general question 
of classifying all such systems is still largely open. In particular, major porgress would be achieved in 
our understanding of the subject if we could obtain a classification or a characterization of all families of 
SL-OPS. (Recall the classification performed by Bochner [2] and Lesky [21] deals only with the classical 
OPS.) It seems clear by now that the Darboux transformation will be one of the key necessary tools in 
achieving such a goal. It should be noted that when referring to the Darboux transformation, we do not 
mean here the factorization of Jacobi matrices into upper triangular and lower triangular matrices |20| . 
Indeed, while such a transformation is defined for any OPS, the transformed OPS will in general not be 
a SL-OPS even if the original OPS was one. We will rather use algebraic Darboux transformationsij, also 
known as rational factorizations, which are defined only for SL-OPS. In these transformations, it is the 
second order operator T that needs to be factorized as the product of two first order operators T = AB, 
and the transformed operator T is obtained by reversing the order of the factors, T — BA. We shall 
see that by construction, these algebraic Darboux transformations transform an SL-OPS into another 
SL-OPS. 

We are now ready to state the main result of our paper: 

Theorem 1.1. Every y^m orthogonal polynomial system for m < 2 can be obtained by applying a sequence 
of at most m Darboux transformations to a classical orthogonal polynomial system. 

The proof of this theorem is done in several steps. The first step, carried out in Section [3l consists in 
the classification of X2 flags and the determination of the corresponding pole structure for the coefficients 
of the second order linear differential operators that preserve them. This forms the substance of Theorem 
13.21 It should be noted that in contrast with the codimension one case, the canonical codimension two 
flags contain free parameters (flag moduli). In Section 5 we provide the necessary background to select 
from the classification of X2-flags those that give rise to a well defined SL-OPS. This selection is performed 
in Section 6, where Theorem 16.11 provides the classification of X2 orthogonal polynomial systems. It is 
worth noting that this classification contains new families of X2-Laguerre and X2- Jacobi polynomials; for 
example the new Laguerre-type family of Section 6.5.6 with weight e~^x^/''/(4x + 3)^. The second step 
in the proof of Theorem 11.11 , which is carried out in Section 21 consists of the proof of the key property, 
stated in Theorem 14. 2[ that all Xi and X2 operators are related to a classical operator by a Darboux 
transformation or a sequence of two Darboux transformations. 

Finally, we will conclude by stating our general, yet-to-be proved, conjecture, which extentcnds the 
result of Theorem 11.11 to arbitrary codimension. 

Conjecture 1.1. Every X„i orthogonal polynomial system for any codimension m can be obtained by 
applying a sequence of at most m Darboux transformations to a classical OPS. 

2. Preliminaries and definitions 

2.1. Polynomial flags. Let V denote the infinite-dimensional space of real, univariate polynomials, and 
let Vn C 7^ be the n + 1 dimensional subspace of polynomials having degree n or less. We define the 
degree of a finite dimensional polynomial subspace U G V to he 

deg U = maxjdegp : p G 14}. (1) 

Deflnition 2.1. A polynomial flag is an infinite sequence of polynomial subspaces C/i C C/2 C . . ., nested 
by inclusion, such that dimC/^ — fc, and such that degL'^fc < degC/fe+i for all k. The total space of a 
polynomial flag is the infinite-dimensional polynomial subspace 

00 
U=\JUk. (2) 

fe=i 
Deflnition 2.2. Let U C V he an infinite dimensional polynomial subspace. A degree-regular basis oiU 
is a sequence of polynomials {pfc}fe^i such that degpk < degpk+i and such that U = spanjpfc}. 

Proposition 2.1. Let Ui C U2 C ■ ■ ■ be a polynomial flag, U the total space, and {pk}'kLi '^ degree regular 
basis. Then, for all k = 1,2, ..., we have 

Uk = span{pi,...,pfe}. 



A wider class of these transformations has been extensively used in Quantum Mechanics to generate new exactly solvable 
problems from known ones. The subclass of interest to us in the context of OPS consists of the set of transformations that 
preserve the polynomial character of the eigenfunctions. This particular class of Darboux transformations was characterized 
in[S[lO]. 
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Proposition 2.2. Let U d V be an infinite dimensional polynomial subspace. Let Uk d U be the unique 
k-dimensional subspace having minimal degree. Then Ui d U2 d ■ ■ ■ is a polynomial flag whose total 
space is U . 

Proposition 2.3. Let Ui d U2 d . . . be a polynomial flag andlA the corresponding total space. Let Uk be 
as above. Then, Uk = Uk. 

The above propositions show that there is a natural bijection between the set of polynomial flags and the 
set of infinite-dimensional polynomial subspaces. Going forward it will sometimes be more conveninient 
to specify the total space rather than the actual flag. The identification of the flag and its total space 
will be implicitly assumed. We will use the complete notation U : Ui d U2 d . . . for the flag and its total 
space, but we will write only U to denote the flag Ui d U2 d . . . where no confusion can arise. 

Definition 2.3. Given a polynomial flag U : Ui d U2 d . . . , we define the degree sequence {f^fc}^l and 
the codimension sequence {w.fc}^-^ as 

Uk = deg Uk, nik =nk + 1 - k. (3) 

where mk is the codimension of Uk in Vn^. ■ 

It is easy to see that {uk} is strictly increasing while {mk} is non-decreasing. In this paper we will 
focus on fiags with finite codimension, which means that the total space U has finite codimension in V, 
or equivalently, that the codimension sequence {ruk} admits an upper bound m — maxkmk, which we 
call the codimension of the flag. As mentioned in the Introduction, one might also characterize m as 
the number of gaps in the degree sequence. We will say that a polynomial fiag has stable codimension if 
wLk = m for all fc, or equivalently if the degree sequence satisfies ni — m and Uk+i = Uk + 1 for all fc > 1. 

The simplest of all polynomial flags in the standard flag U^t '■ Vq d Vi d 1^2 C . . . . The total space for 
this fiag is V, its degree sequence is NU {0} and it has stable codimension zero. 

Definition 2.4. We will say that a second order differential operator 

T[y(z)] = p{z)y" + q{z)y' + r{z)y, (4) 

is rational, if the coefficients p, q, r are rational functions of the independent variable z and the prime 
denotes derivation with respect to this variable, y' — -^- The poles of a rational operator T are the poles 
of p, q and r. An operator T with no poles is said to be polynomial. If there is one or more poles, then 
we will refer to T as non-polynomial. 

Definition 2.5. We say that a polynomial flag U : Ui d U2 d . . . is invariant under a rational operator 
T[y] if T{Uk) d Uk for all k. We let T>2{hl) denote the vector space of all second order operators that 
preserve the flag Uu 

In the analysis of invariant polynomial flags, no generality is lost by considering only second order 
operators with rational coefficients, as evidenced by the following 

Proposition 2.4. Let T[y] = py" + qy' + ry be a differential operator such that 

T[yi\ ^ gt, i = 1,2,3, 

where yi,gi are polynomials with 2/1,2/2,2/3 linearly independent. Then, p,q,r are rational functions. 

Proof. It suffices to apply Cramer's rule to solve the linear system 

'y'{ y'l vi 

y'i y'2 V2 

\y'i y'3 ys 

D 

Definition 2.6. A polynomial flag is imprimitive if it admits a non-trivial common factor. Otherwise, 
the flag is said to be primitive. 




We stress that invariance of the whole flag W : C/i C f72 C . . . is a much stronger condition than the invariance of 
the total space ti. For the purpose of this study, we will always require invariance of the flag, since this condition leads to 
polynomial eigenfunctions of the operator. 
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Proposition 2.5. Let U be a primitive flag, let ^ be a polynomial of degree > 1 and let 

U = ^jU = {fip -.peU}. 

be the corresponding imprimitive flag. Suppose that T[y] is a rational operator that preserves U. Then, 
the gauge-equivalent rational operator T — ^T^^^ preserves lA . 

Therefore, primitive flags can be regarded as canonical representatives for the equivalence relation modulo 
gauge transformations, and we can restrict our attention to primitive flags in the classification of invariant 
polynomial flags. The main object of our study is then the class defined in the following definition. 

Definition 2.7. A second order operator that preserves a primitive polynomial flag, but does not preserve 
the standard flag will be called an exceptional operator. An exceptional flag is the maximal primitive 
polynomial flag that is preserved by a second order exceptional operator. Exceptional flags and operators 
of finite codimension m> 1 will henceforth be called X^ fiags and operators. By contrast, a second order 
differential operator that preserves the standard flag V, will be referred to as a classical operator. 

Theorem 2.1 (Bochner). A classical operator has the form 

T[y] = py" + qy' + ry 

where p (^ V2, q (z Vi are polynomials of the indicated degree, and where r is a constant. 

Proposition 2.6. An exceptional operator is, necessarily, non-polynomial. 

Proof. See the proof of Lemma 3.1 in [16]. D 

Thus, an exceptional operator has poles, but it also has an infinite number of polynomial eigenfunctions. 
When classifying exceptional fiags by increasing codimension, each flag will give rise to new operators not 
considered at lower codimension, which justifies the definition above. Here are some examples to illustrate 
these definitions 

Example 2.1. The flag with basis {1, z^, z^, . . . , } is exceptional because the operator 

T[y]^y"~^-^ 

preserves the flag. The degree sequence is {0, 2,3,...} and the codimension sequence is {0, 1,1,...} so 
the fiag has non-stable codimension 1. 

Example 2.2. By contrast, the fiag spanned by {z + 1, z^, z^, . . .} has a stable codimension m — 1. This 
flag is exceptional because it is preserved by the operator 

Example 2.3. Let Hk{z) denote the degree k Hcrmite polynomial. The codimension 1 flag spanned by 
{Hi, H2, H3, . . .} is not exceptional. The fiag is preserved by the operator T[y] — y" — zy' . However, this 
operator also preserves the standard fiag, which violates the maximality assumption. 

Example 2.4. The codimension 1 flag spanned by z, z^, z'^, . . . is preserved by the operator 

T[y\ =y + ^. 

Z Z"^ 

This is not an exceptional flag because it is imprimitive as z is a non-trivial common factor. In fact, the 
operator T is gauge equivalent T = zTz~^ to the operator T[y] — y" that preserves the standard flag. 

Example 2.5. Let 

y2fe-i=22fc-i-(2A:-l)z 

J/2fc = z^'' - kz 

Consider the fiag spanned by {1, 1/3, 2/4, J/s, . • ■}■ The degree sequence of the fiag is 0, 3, 4, 5, . . . so it is a 
non-stable codimension 2 flag. The fiag is preserved by the following operators : 

Ts[y] = {z^^l)y"~2zy', (6) 

T2[y]^zy"-2^l + -^^y\ (7) 

Ti[y]=y" + z(l--^)y'. (8) 



2k U.2 fc = 2,3,4,... (5) 
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The flag is exceptional, because Ti and T2 do not preserve the standard flag. Since T2,Ti have 2 distinct 
poles, they do not preserve a codiniension 1 flag (see Lemma l3.3p . 

3. Classification of exceptional codimension 2 polynomial flags 

In this Section we perform a classiflcation of all X2-flags up to affine transformations of the independent 
variable z. We exhibit degree- regular bases for each of them, and we determine the X2-operators that 
preserve them. We begin by introducing the following flags: 

£'^'\a;b)~{per:p'{b) = ap(b)} (9) 

f'"'(ao,ai;6o,6i)— £*''(ao;6o)n£:<''(ai;6i) (10) 

£*^'(aoi,ao3,a23;&) := {p £ P : p'{b) = aoip(b), p'" (b) = 3a23p"(6) + 6ao3p{b)} (11) 

The first fiag has codimension one and its associated Xi-operator will have a simple pole at z = &. The 
second flag has codimension two and its associated X2 operator will have two simple poles at 60 and 61. 
The third flag has codimension two and its associated X2 operator will have a simple pole at b. The 
notation in the superindices is connected to the the order of the poles of the weight for the exceptional 
orthogonal polynomial system based on the flag. This will become clear in Section [6l In any case, the 
sum of superindices must always coincide with the codimension of the flag. 

Some, but not all of the parameters in the above flags can be normalized by means of an affine 
transformation. Thus, unlike the codimension one case, the X2 flags contain free continuous parameters, 
which shall be refereed to as flag moduli. As explained before, the parameters 6, 60 and 61 will be the 
positions of the poles of the operators. If there is one pole we will set 6 = and if there are two poles we 
will normalize them as 60 = 0, bi — 1. Note that any two poles in the complex plane can be transformed 
into and 1 by a complex affine transformation, so there is no loss of generality involved in the above 
normalization. 

Below, we describe each of the above flags in terms of a basis. 

£:'^'(a;0) =span{l + az,^^^^^■*,...} (12) 

£:<">(ao, ffli; 0, 1) = span{2''((ai - 2){z - 1) + 1), (z - if^ao + 2)z + 1)}U (13) 

{2^(2 - l)^2^}°lo, 

£*^'(aoi,ao3,ffl23;0) = span{l + aoiz + ao3Z'^,z^ + a23Z^ , z"' , z^ , . . .} (14) 

Let us first recall the main result of the classification of Xi-flags first proved in [12] (see [16] for a more 
recent and streamlined proof). 

Theorem 3.1. Every stable Xi polynomial flag is affine- equivalent to 

f (1) (1; 0) = span{l + z, z^ z^ , z^ . . .}. 

Every unstable Xi polynomial flag is affine- equivalent to the monomial flag 

f(i)(0;0) =span{l,z^^^z^...}. 

Note that, as mentioned before, the most general Xi flag up to affine transformations contains no flag 
moduli. The main result of this section is the following theorem that describes the situation for X2 flags. 

Theorem 3.2. Up to an affine transformation every X2 flag is equivalent to one of the following two 
flags: 

(1) £(")(ao,ai;0,l) 

(2) £(^'(001, ao3, a23; 0) subject to the constraint 

ao3(aoi - a23)(6ao3 + aoia23(aoi + 023)) = (15) 

Before we can address the proof of the above theorem, we need to introduce more concepts and establish 
some key intermediate results. 

For a polynomial y{z) and a constant 6 G C, we define ord^ y > to be the order of 6 as a zero of y{z). 
Let U C V he a, polynomial subspace. For 6 e C define 

Ib(U) = {oTdby:yeU}. (16) 
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Lemma 3.1. Let T be a rational operator that preserves a primitive polynomial subspace U dV. Let 

oo 

T= j2t,, 

i— — d 

where 

Ti[y] = z* {piZ^y" + q^zy' + r^y) 

for some constants Pi,qi,ri be the degree-homogeneous representation of T in terms of Laurent series. If 
T has a pole at z = 0, then d = 2, r_2 = 0, and there exists a positive integer a > 1 such that 

/o(W)=N/{l,3,...,2a-l} (17) 

Proof. Observe that Ti is degree-homogeneous, meaning that 

Uz']^{p^J{J-l)+q^J+n)z'+^. 

So either Ti annihilates a given monomial z^ , or it shifts its degree by i. A non-zero Ti can annihilate at 
most two distinct monomials, whose exponents j satisfy the quadratic constraint 

Ptjij - 1) + qi.1 + ri =0. 

By definition, i G /o if and only if the flag contains a polynomial of the form z*+ higher degree terms. 
Since T preserves U and since T_d is the leading term of the operator, it follows that T_d preserves the 
monomial subspace {z^ : i £ Iq}. 

For T to have a pole at z = we must have d > 0. Since U is primitive, G /o and therefore T-d 
must annihilate z° = 1. Observe that the leading order d must also be d > 2, since d = 1 would require 
that r_i[l] — ^ r_i — 0, so operator T would be polynomial, contrary to the hypothesis. To conclude 
the proof, we will establish that d has to be precisely 2. Since the fiag U has finite codimension, there 
are only a finite number of gaps (missing integers) in the set Iq . Let i ^ Iq he one such gap, then either 
i + d ^ Iq, or T_ci annihilates z*+'^. Hence, 1 ^ /o must be a gap. Otherwise, since d > 2, T_d would 
need to annihilate three monomials: z°, z^ and at least one higher degree monomial, which is impossible. 
Thus, for some integer a > 1, the gaps in the /q sequence are 1, 1 + d, 1 + 2d, . . . , 1 -I- d{a — 1) ^ /qj with 
T-d[z'^°'^^] = 0. Note that T_d annihilates 1 and z'^°'~^^ so it cannot annihilate any other monomial and 
therefore the above gaps are the only possible gaps in Iq. It follows that 2 G /q is not a gap. If the leading 
order was d > 2 then T-d would be required to annihilate also z"^, which is impossible. We conclude then 
that d = 2 and since r_2[l] = we must have r_2 = 0. The assertions of the lemma are proved. D 

The following lemma shows how to decompose a rational second order operator that preserves a prim- 
itive polynomial fiag. 

Lemma 3.2. Let T be a second order rational operator with poles 6i, . . . , fe^r € C. If T preserves a 
primitive polynomial flag of finite codimension, then necessarily it has the form 

N 

T[y] = p-2y" + {p-izy" + q-iy') + {paz^y" + qoy' + ray) + ^ Q ■ 



y' - a,y 



z - bi 



where pi,qi,ri G K. and a^, q G C are constants. 



We see therefore that an exceptional operator must have rational coefficients that can only contain 
simple poles. 

Proof. We decompose the given operator as 

N 

T = y^ r*^*^ 

i=0 

where T*-^-* is a polynomial operator and where 

rWf 1 = '^-i^ q-kz - h)y' + r^:!ly ^ Pijjz - bify" + qijjz - bi)y' + r^^y 

^y^ z-bi iz-b,y +^4^ {z-b,y 

for some positive integer di > 1 and constants Pij , qij , rij . 
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Let U be the total space of the primitive flag preserved by T. Since T{U) C V, it fohows that 
T^'^'^{U) C V for every i = 0,1, ... ,N. By construction, the operators T'^^\ . . . , T'-^^ ah fower degrees. 
Since T preserves an infinite flag, it cannot have a degree raising part. Therefore, T^*^' has the form 

T(°)[y] = p_22/" + ip^izy" + q-W) + (po^'y" + go2/' + r^y) 

Expanding the operator coefficients as Laurent series in z — 6i, we apply Lemma l3. II to conclude that 
di = 2, r_2 = for alH = 1 . . . , A^. The desired conclusion has been established. D 

Note that if bi is a real pole, then the constants a^ and Ci must also be real since the flag is real too. 
The next lemma shows that for every pole bi of an exceptional differential operator, the elements of its 
invariant flag must satisfy a first order differential constraint at that pole. 

Lemma 3.3. Let T[y\ he a second order rational operator with poles bi, . . . ,1)^ that preserves a primitive 
flag lA of finite codimension. Then, there exist constants ai, . . . , an such that the elements of y Cz U obey 
1st order differential constraints of the form 

y'{bi)^a^y{bi), i^l,...,N. 

Proof. By Lemma [O] for each i = 1, . . . ,N the total space U contains a polynomial of the form 

y^o'\z) = l + a,{z-b,) + 0{{z-b,f), 

but does not contain an element of the form 

{z-b,)+0{{z-b,f). 

Therefore, every y ^ U either starts as y^ (z) or at degree 2 in (2; — bi) so in any case it obeys the 
constraint y' {hi) — aiy{hi). D 

At this point, it becomes necessary to describe and analyze certain degenerate subclasses of the S^^^'' 
and £^^^ flags defined in (fTO |) -(fTT |) . The distinguishing property of these subclasses is the first two elements 
of the degree sequence of the flag. Thus, when we write £ij where < i < j < 3 we are referring to a 
codimension two flag whose degree sequence is {i, j, 4, 5, 6, ... }. The generic case is the stable codimension 
two flag £23, which starts at degree 2 and has polynomials of all degrees fc > 2. We analyze each of the 
above 3 families in more detail, and then give a proof of Theorem 13.21 

Proposition 3.1. The £^^^' flags are classified into the following subclasses, according to their degree 
sequence: 

^23 = ^ (flOi cti; 0, 1)1 'with fliflo + fli — flo 7^ (18a) 

= span{(aoai + ai ~ ao)z'^ + (2 — ai){aoz + 1), {z — 1)^(1 + (2 + ao)z)} U 
span{z2(z _ l)2zJ'}=^„ 

^1"^ = ^*"^ f «o, T^^; 0, 1) , wtth ao + -1, and (ao, a^) i {(0, 0), (-2, 2)} (18b) 

\ t + ao / 

= spanjaoz + 1, (2 - 1)^(1 + (2 + ao)z)} U span{z2(z - \-fz^}f=Q 
f^Ji) =f(")(0,0;0,l)==span{l,(z- 1)2(1 + 2z)} U span{z2(z - l)2z^}°^o (18c) 

f (11) = f (11) (-2, 2; 0, 1) = span{2z - 1, z^} U span{z2(z _ \fz^}f^^ (18d) 

Proof. This follows by direct inspection of ([T3|. D 



Also note that £\^2 can be obtained as a limit of f 23 by setting oq — t — 2, ai — t + 2 and then sending 



i = 0. The flags in Proposition 13.11 are all X2-flags whose operators have two simple poles at and 1. In 
the following Proposition we provide a basis for the I?2-spaces of operators that preserve them. 

Proposition 3.2. The generic flag £33 'las a 2-dimensional T>2 space. The non-stable flag f^g has a 
3-dimensional T>2, while £^^ and £^2 both have a 4-dimensional T>2. The most general second order 
operator that preserves each of these flags is shown below (and therefore a basis of their T>2-space). The 
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symbols ao,ai denote the flag moduli while the symbols c,Co,ci,qQ, X denote free constants appearing in 
the operator 



T, 



(11) r 



-z (ao — ai)(ao — ai + 4) — z{a^)al — ao — fli + 3ai) 



al 



ai\y + 



n{ll)r 



+ c (2((ao - ai)(aoai - 2ao + 2) + 2ao) + (ao - l)a\ - (ao - 3)ai + ao(ao + 1)) j/' + 

cao(ao + 2) , c(ai - 2)ai , 
H ^ _ -^ [y - aiy) + [y - aoy) + Xy 

= f -(cQ + ci)— + Co ( 2 - -^ ) ) y" + ((aici - aoco)2 + (ao - l)co + ci)y' + 

— {y'-aoy)A \{y' - aiy) + Xy, ai = ^ 

z 2: — 1 ao + 1 



(19a) 



(19b) 



n(ll)r 



n(ll)r 



I / \ z qoZ ( 1 

= ( - (go + CO + Cl) y + — + CO I 2 - - 



y" + ( go ( 2 - - ) - Co + ci ) J/+ (19c) 



Cl 
2-1 



y +Xy 



= ( (co + Cl - go) Y + (^y - cij 2 - y j y" + f go f 2 - - ) - 2co + 2ci ) y'+ 



+ ^{y' + 2y) + ^{y'-2y) + Xy 

2 2—1 



(19d) 



Before we turn to the proof of this last Proposition, observe the duahty between flag moduli and free 
parameters in the operator. In the general case £"23 the flag has two moduli (ao, ai) and the I?2-space has 
dimension two. In the case £"^3 the flag has one modulus ao and the operator has three free parameters, 
since diml?2 [^13 ) — 3. In the last two cases £q^ and £{2 the flag is completely specified (no flag 
moduli) but the operator contains four free parameters. 



Proof. By Lemma 13. 2[ we must consider an operator of the form 

rpr 1 // , / // , /N , . 2 // , /N , co{y'-aoy) ci{y' - aiy) 

T[y\:^p~2V +{p-izy + q-iy ) + {poz y + qoy ) -\ h , 

z z — I 

where p_27P-i7'Z-iiPoi 901^1,00 are undetermined coefficients that need to be constrained so that T 

preserves the flag in question. Applying the relation 

y'(0) = ao2;(0), yeU 

to the constraint 

r[y]'(0)-aoTM(0) = 0, yeU (20) 

yields: 

(co/2 + p_2) ?/"(0) + (3aoP-2 - (3ao/2)co - ci + p_i + q^i) y"{0)+ 
+ (floCo + (ao - ao + ai)ci - ao^-i + ao^o) y(0) = 

Since y"'{0),y"{Q),y{0) vary freely for y eU the coefficients of all 3 terms must vanish in order for (PH)) 
to hold. An analogous constraint holds for 

T[y]'{l)-a,T[y]{l)^0. 

Since there are 7 parameters and only 6 linear, homogeneous constraints, there exists at least one non- 
trivial operator that preserves S'^^y The desired solution vector 

[P-2,P-i, q^i,Po, qo, Cl, coY 



belongs to the null-space of the following matrix 

/ 1 

-ao 1 1 







\ 



(21) 



1/2 

-3ao/2 

al 



1 

\0 -af -(ai - l)ai af ag- ai{ai + l)J 

A direct calculation shows that all 6 minors of the above 6x7 have aoai -|- ai — ao as a factor, and that it 
is not possible for all the minors to vanish if aoai -|- ai — ao 7^ 0. Hence, generically the above constraint 



-af, 

110 1 

-ai 1 — ai 1 2 — fli 



do 

1 



-1 

al- ao + ai 
1/2 

-3ai/2 

„3 
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matrix has rank 6, and there exists a unique, up to a scalar factor, solution, which after some calculation 
provides the operator Tjg . 

Setting oi = ao/(ao + 1) in the above matrix drops the rank of the matrix to 5, provided, Oq ^ {0, —2}. 
Now the nullspace is 2-dimensional; this gives the form of T-^-g . Setting oq = ai = in the constraint 
matrix gives a matrix of rank 4. The nullspace corresponds to the operator Tqj . Similarly, oq = — 2, ai = 
2 also gives a rank 4 matrix, whose nullspace corresponds to the operator T^^ ■ ^ 

Proposition 3.3. The flag E^^ is an X2 flag, provided Uq ^ {0,-2} and ai ^ {0,2}. The non-stable 
flags £[l^\£f^l^\£[l^' are all X2 flags. 

Proof. It is clear that all the operators preserve codimension two flags and since they have poles they 
do not preserve the standard flag. It only remains to prove the maximality assumption, i.e. that these 
operators do not preserve a flag of codimension one. By Lemma |3.3[ an operator with two poles cannot 
preserve a codimension 1 flag. By inspection, if ao,ai satisfy the conditions given above, the operator 
T23 cannot preserve a codimension 1 flag. On the contrary, if oq = 0, a direct calculation shows that 
^23 [1] =0 and hence £"'^'(0,01) is not the maximal flag preserved by T23 . Similarly, if oq = —2 then 
2z — I is again an eigenpolynomial. Similar remarks hold for the cases ai = and ai = 2. 

For the degenerate, non-stable flags, by taking cq, ci 7^ we obtain operators that preserve these flags, 
but have 2 distinct poles. Therefore, by the same argument, these operators cannot preserve a flag of 
smaller codimension. D 

We now turn to an analysis of the one-pole X2 flag f'^-* . In the language of Lemma [5TT1 this flag is the 
most general codimension 2 flag with the order sequence /o = {0, 2, 4, 5,6,.. .}. The Lemma below derives 
the constraint ([T5|) as the necessary and sufficient condition for such a flag to have a non-trivial T>2. 

Lemma 3.4. Every X2 flag that is preserved by an operator with a unique pole is translation-equivalent, 
to £'^^'(aoi, ao3, 023; 0) where the parameters satisfy (jI5p . Up to a multiplicative constant, a second order 
operator that preserves .such a flag has the form 

T(2) [y] = y" + {p.,zy" + q.,y') + [p^z^' + q,y') - 4 ^^''"°^^^ + Ay (22) 

z 

where 

V-i = 2aoi - 2a23 (23) 

g_i = -7aoi -I- 5a23 (24) 

Sagi - 6ao3 - 5aoia23 
ao3(aoi -a23)(aoi +023) | | 9o | = ( ) (25) 

6ao3 + 5aoia23 - 3033 

Proof. By Lemma 13.21 an operator with a unique pole at z = that preserves a polynomial flag has the 
form 

T[y] = p^,y" + {p^.zy" + q^,y') + {p.z'y" + q^y') + M^lJ^ + Ay 

where p-2,P-i, (Z-ijPo, 90, A, c are undetermined coefficients. If we demand that the flag has codimension 
2, then the flag must be £(^^. By Lemma [5TT1 it follows that we must also require T_2[2;^] = 0, where 

T-2\y\ ^p-2y" + — , 

z 

This imposes the condition c = — 4p_2 and since we require a non-trivial T_2, we must have p_2 7^ 0- 
Hence, without loss of generality, we impose 

c = -4, p-2 = l 

from here on. The flag £^ in pT|) is defined by the first and third order conditions 

y'(0) = aoiy(O) (26) 

y"'(0) = 6ao3y(0) + 3a23y"(0) (27) 



and where po,qo satisfy: 






a-oi 


2ao3 


ao3 


4a23 


0.23 
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Imposing these conditions on T[y] yields 

(5aoi - 3a23 +P-i +q-i)y"{0) + (-4aoi - 6ao3 - aliQ-i + aoigo) ?/(0) = 

(floi + 023 + 3p_i + q-i) y^^H0) + 

+ 3 (-4aoia23 + 6ao3 - 6033^-1 - ial^q^i + 4a23Po + a23go) y"{0)+ 

- 6ao3 (4aoi + 4aoia23 + 6a23P-i + (3a23 + aoi)?-! - 6pa - iqo) y(0) = 

The values of j/'-'*'(0), j/'-^'(0), j/(0) vary freely for y <E £'--^\ and hence, invariance holds if and only if the 
coefhcient of each of these expressions vanish. The conditions ([23|) (p4|) follow from the vanishing of the 
leading order coefficients. Once these values of p_i,g_i are imposed, the overdetermined constraint (P5|) 
expresses the vanishing of all the remaining coefficients. The vanishing of the determinant of the matrix 
in (|25p is the compatibility condition for these constraints, and this is precisely condition ([T5|) . D 



As we did before for the two-poles X2 flags, the one-pole X2 flags can be classifled according to their 
degree sequence. 

Proposition 3.4. Every one-pole X2 flag is affine- equivalent to one of the following: 

f}f^(a) :=£(2)(l,0,a;0) = span{l + ^,^2^az^^^z^...}, a^O, (28a) 

Sg"'^ :^ f (2)(o, 0, 1; 0) = span{l, z^ + zJ, z^ z^ . . .}, (28b) 

^{2a) ,_ £(2)(i^ 0, 0; 0) = spanjl + z, z^, z^ z^ z^ . . .}, (28c) 

£:^2a) .^ £(2)(Q^ 0, 0; 0) = span{l, z^, z^ z^ z^ . . .}, (28d) 

£if\a) ■.^£'-^\a,a,a;0) = spa.n{l + az - z^ , z^l,! + az), z^ , z^ ,.. .}, a ^ 0, (28e) 

£^f\a):=£^^'>{a,-a{a + l)/6,l), a^Q (28f) 

= span{l + az + a{a + l)z^/6, z^ + z^, z"*, z^, . . .} 

Proof. The three types of flags labelled (2a), (25) and (2c) correspond to the three different ways of 
satisfying the defining constraint (fT5|) on the three flag moduli. The type (2a) flags are obtained by 
applying the constraint ao3 = to a general type f'-'^' flag. By (I14p . the resulting degree regular basis is 

1 -I- aoiz, z^ + a23Z^, z'', z^, . . . 

If aoi,a23 7^ 0, then a scaling transformation can be used to send one (but not both) of the above 
parameters to 1. The various subclasses listed above arise if one or both of aoi, a23 — 0. 

The type (26) flag is obtained by applying the constraint a23 = ffloi- An examination of (|25p shows that 
it is not possible for a23 = aoi — 0, ao3 7^ 0. Therefore, for the type (26) subcase, we must have aoi 7^ 0. 
Thus, in this case, transforming (J14p to a degree regular basis gives 

1 I „ y "03 2 ^2 I ^ ^3 ^4 ^5 

1 -f- agiz — z , z -T aoiz , z , z , . . . 

aoi 

Now a scaling transformation can be used to set ao3/aoi — 1. 
The type (2c) flags are obtained by imposing 

am = -aoia23(aoi + 023)76. 
In this case, the degree regular basis is 

1 + aoiz + aoi (aoi + a23)z^/6, z^(l + a2'iz) , z^ , z^ , . . . 

No generality is lost if we assume that aoi, a23 7^ 0, because otherwise we will obtain a flag of type (2a). 
Finally, a scaling transformation is used to set a23 = 1. D 

Note: as above the flag subscript indicates the degree sequence of the flag. 

Proposition 3.5. The flags ^13° , 1^23 '^23"^ have a 2- dimensional 2?2- The degenerate flags £q^\£\^^ 
have a 3- dimensional V2, while Eq-^ has a 4- dimensional Vi- The most general second order operator 
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that preserves each of these flags is shown below. The symbol a represents the flag modulus while the 
symbols c,pQ,qo,X are free constants that appear in the operator. 



= c m - 3a) (3 - a) ^ + 2 (1 - a) z + 1 j y" + 



t[T^ M = c ( (1 - 3a) (3 - a) ^ + 2 (1 - a) z + 1 ) y"+ (29a) 

+ c ((5a - 3) a^ + 5a - 7) y' + ^^^^ ~ ^ ' + Ay 



-,(2a)r 



('(3c - go) \ + c(l - 2z)\ y" + (5c + qoz) y' 



= ( (3c - go) ^ + c(l - 22) ) y" + (5c + qoz) y' - ^ + Xy (29b) 

Ti'^"' [y] = (P02' + c{2z + 1)) y" - c (32 + 7) y' + ^S^l-jll + Ay (29c) 

T!n\y] = (P"^' + c) y"(2) + go2y'(2) - ^^^ + Ay (29d) 

TS'\y] = c[l^z^ (a^ + 3)) y" + c (22 (a^ + 3) - 2a) y' {z) + ^''^"^^~ ^'^ + Ay (29e) 

rir'M = c(l + (a - l)2)2y" + c{{a - 1)(1 - 3a)2 + 5 - 7a)y' + ^'"'^'^^ ^ ^ ^ + Ay (29f) 

Proof. Each of the flags in question is a speciahzation of the f^^^ flag discsussed in Lemma [3.41 imposed 
in such a way so that (fTSl) holds. The 3 factors in dTS]) give us the 3 possible cases: £^'^°'\ £^'^^\ E^'^'^\ 
Imposing the respective constraints 

ao3 = 0, a23 = aoi, ao3 = -aoia23(aoi + a23)/6 

transforms (|25p into a consistent, rank 2 system. We can further ehminate one more parameter by means 
of an appropriate scahng transformation. The form of the operators shown above follows from (|23p ([M)) 
and the solution of the corresponding (P5|) . D 

Proposition 3.6. The flags ^13" , i^oa" , 1^12° j '^02° ' '^23 '^23^^ '^'"'^ '^^^ -^2 fl^igs. 

Proof. For each of the above flags, we have exhibited a singular operator that preserves it. It remains 
to show that these operators cannot preserve a flag of smaller codimension. By Lemma 13.11 an Xi flag 
preserved by an operator with a pole at z = 0, must have elements of order 0, 2, 3,4,.. .. Therefore, it 
suffices to check that T_2 (see the Lemma for the explanation of the notation) does not annihilate z^. For 
each of the operators shown in the preceding Proposition, 

T-2[y\ = y • 

z 



Hence, 



r_2f2^1 = -6z. 



-2[ 

Therefore, none of these operators can preserve an Xi flag. D 



Proof of Theorem \3.2l By the above Lemmas, an X2 operator has either one or two poles. In the last 
case, the corresponding X2 flag satisfies two distinct first order condtions 

y\bi) = aiy{bi), i = l,2 

Applying an affine transformation, no generality is lost if we assume that the poles are at 2: = and z ~ 1. 
This gives us flags of type £*^^^). The corresponding X2 operators are given in Proposition 13.21 The X2 
assertion is verified in Proposition 13.31 

In the case of one pole, without loss of generality the pole is at z = 0. In this case, the flag satisfies a 
first and a third order condition, which gives us a flag of type f "'. As it was shown in Lemma 13.41 the 
moduli of the general flag must satisfy the constraint p^ . This gives us the three cases: £^(2") , f (^^) , ft^"^) . 
The corresponding operators for these flags are given in Proposition 13.51 and the X2 condition is verified 
in Proposition l3.6l D 

4. Factorization of exceptional operators 

The results in this section are concerned with factorizations of the differential operators that preserve 
X2 flags and their connection to the Darboux transformation. The usual Darboux transformation involves 
Schrodinger operators and square-integrable eigenfunctions but for our purposes it will be convenient to 
generalize it to second order operators with rational coefficients. 
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Definition 4.1. Let T be a second order differential operator that preserves a polynomial flag U. Let 

T = BA + Xo (30) 

be a factorization of T where A, B are first order operators with rational coefficients and Aq is a constant. 
If the partner operator defined by 

f = AB + \fi. (31) 

also preserves a polynomial flag U we will say that T and T are related by an algebraic Darboux transfor- 
mation. 

Definition 4.2. More generaly, we will say that two operators T and T are Darboux- connected if there 
exists a sequence of algebraic Darboux transformations that connect them. 

The same notion can be defined for polynomial fiags in the following manner: 

Definition 4.3. Two polynomial fiags U : Ui C U2 C . . . and U : Ui C U2 C . . . are Darboux- connected 
if there exists two first order rational operators A and B such that one of the following three possibilities 
occur: 

A[U,] C U^, B[U,] c f/„ t> 1; (32) 

A[U,+i](lU,, B[U,]ciU,+i, i>l, A[Ui]=0; (33) 

B[U,+l]cU^, A[U,]cU,+i, i>l, B[Ui]=0. (34) 

In accordance with [M] we will refer to the above cases as formally isospectral, formally state-deleting 
and formally state-adding. 

Note that this implies that the second order operators T — BA and T = AB preserve the fiags U and 
U respectively, so Darboux-connected polynomial flags are always invariant. It is common to refer to the 
operators A, B as intertwining operators, or simply as intertwiners. 

Definition 4.4. We say that a polynomial flag U is an m-step flag if there exists a sequence of m Darboux 
transformations that connect U to the standard flag. 

Our main results in this section are summarized in the following two theorems: 

Theorem 4.1. Every Xi flag is a 1-step flag. Every X2 flag is either a 1-step or a 2-step flag. 

Theorem 4.2. Every Xi and X2 operator is Darboux-connected to a classical operator. Furthermore, the 
intertwining operators that connect the classical operator to the ^-operator also connect the standard flag 
to the exceptional flag. 

As we show in the next section, one consequence of Theorem 14.21 is that all X2 and Xi orthogonal 
polynomials can be expressed as certain Wronskians involving classical OPs. 

Using the classification of Xi and X2 flags from the preceding section, the proof of Theorem 14.11 is 
broken up into a series of Lemmas. It turns out that Theorem 14.21 is a consequence of Theorem 14. II Our 
proof strategy is to show that if two polynomial flags are Darboux-connected, then so are the operators 
that preserve them. This fact is established by Lemmas 14.21 14.41 and 14.51 We complete the proof of 
Theorem l4.2l at the end of the present section. 

Lemma 4.1. Every Xi polynomial flag is a 1-step flag. 



Proof. Let U — £(a;b) be an Xi flag as per Theorem 13. II Without loss of generality, 6 = 0. Define the 
1st order operators 

A[y]:=l^^, B[y]:=zy'-{az + l)y (35) 



A[U,]cr,-i, 1 = 1,2,... 
B[ynO) - aB[y]iO) = y'(0) - ay(0) - y'{0) + ay{0) = 



By inspection. 

Also, 

Hence, 

B[V^-l](lU„ i = l,2,.... 
Therefore, AB preserves the standard flag, while BA leaves invariant Ui for every i = 1, 2, . . . . D 

Lemma 4.2. Every Xi operator is Darboux-connected to a classical operator. 
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Proof. Let U — £{a; b) be an Xi flag as per Theorem 13. II Without loss of generahty, 6 = 0. Let 

Aai [y] = A[y] + aiy', B^^ [y] = B[y] + a2Z^y\ 

where A,B are the operators defined in (|35|) . Observe that 

A^,[U^]^V^-l, Bo,,[V^-i]cU,, i = l,2,.... (36) 

and that 

dini{cBQ2^ai + A : «!, a2, c, A e K} = 4. (37) 

It follows that every operator in the vector space in (|37|) preserves the Xi flag. In [IX, Proposition 4.10] it 
was shown that diinl?2(^) = 4. Therefore, by dimensional exhaustion, every operator T € V2(U) admits 
a rational factorization of the form T = cBa2Aai + A. To conclude, we observe that, by p6|) . the partner 
operator T = cAa-^B^^ + A preserves the standard polynomial flag. D 

Lemma 4.3. Let U be a polynomial flag. //diml?2(^) > 2 then there exists a second order operator 
T € 'D2(U). //diml'2(^) = 2, exactly, then {l,r} is a basis of'D2{U). 

Proof. It is clear that 1 G 2?2(^)- If there exists a first order operator S E 2?2(^)i then 5^ G 2?2(^) is a 
second order operator, as was to be shown. It also follows that, if 'D2(U) contains an operator of the 1st 
order, then dimI?2(Z^) > 3. Hence, if dimI?2(W) = 2, exactly, then every T S 'D2iU) is either a constant 
multiplication operator or an operator of the second order. D 

Lemma 4.4. Let hi <Z V be a polynomial flag and A[y] a 1st order operator such that hi :~ A[hi] d V is 
also a polynomial flag. Furthermore, suppose that A[Ui] = {0} and that diml?2(^) > 2. Then, hi, hi are 
Darboux connected. Furthermore, every operator in 2?2(W) is Darboux- connected to an operator in 'D2{hi). 

Proof. Choose a non-zero (j) E Ui. Let T E 'D2{hi) be given. Since </> spans Ui and since T[Ui] C Ui we 
must have 

{T-XM=0 

for some A G R. Write 

T[y] = py" + qy' + ry 
A[y] = b{y' - wy) 

where p{z),q{z),r{z),b{z) are rational functions and where w{z) — (j)' {z) / (j){z) , because A[0] = 0, as per 
the above assumption. Next, set 

B[y] = b{y' - wy) 
where 

w = —w — q/p + b' /b, b = p/b 
A direct calculation then shows that 

T = BA + \. 
Since the kernel of Ajt/i+i is 1-dimensional we actually have 

U,^A[U,+i], 1 = 1,2,... 
Since 

r[t/,]cc/„ 1 = 1,2,... 

it follows that 

B[U,](lU,+i, 1 = 1,2,... 

Therefore AB E V2{hi) and BA E X>2(W). By Lemma |131 there exists a T e X>2(W) such that p{z) ^ 0. 
This proves that hi and hi are Darboux connected. D 

Lemma 4.5. Let hi, hi be Darboux- connected polynomial flags. If dmi'D2{hi) = 2, then every operator in 
'D2{hi) is Darboux- connected to an operator in D2{hi). 

Proof. Let A[y] and B[y] be 1st order operators that connect the two flags. It is clear that T = cBA + A 
preserves hi for all c, A G M. By exhaustion every operator in 'D2{hi) has this form. By assumption, the 
partner operator T = cAB + A preserves the partner flag hi. D 

Lemma 4.6. The flag £23 'is a 1-step flag. 
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Proof. Recall that £33 = ^'■^'' (o-Oi oi! 0, 1) where aoai + ai — ag 7^ 0. Consider the 1st order operators 

z z — 1 

B[y] = z{z - 1)(2 - ai + (ai - oq - 4)z) y'+ 

+ ((aofli + fli - ao)z^ + (2 — ai)ao^ + 2 — fli) y 

Let C/i C C/2 C . . . be the flag corresponding to the total space £33 ; see (|18ap for a degree regular basis. 
A direct calculation shows that 

B[y]'{Q) - aoi?[y](0), B[y]'{l) = a,B[y][l\ 

Since i? raises degree by 2, it follows that 

B[Vj-i]cU,, J = 1,2,... 

From the definition (fTO|) . we see that ^[£33 ] C V. Furthermore, 

((ai - ao)j + aoai)z-'"i - z^~'^jai 



A[z^] 



z- 1 
(oi - j)aQ 



(floai + j(ai — ao))z^ ^ + lower deg. terms 



z-1 
Since deg C/j = j + 1, it follows that 

A[[/,] c T'.-i, j = l,2,... 

as was to be shown. D 

Lemma 4.7. T/ie flag £{^ is a 2-step flag. 

Proof. The degree regular basis is shown in (|18bl) . In particular, 

Ui = spanjl + aoz}. 

Define 

aiW[y,l + aoz] y' - a^y y' - aiy ao 

A[y\-^ n — ^ — ""1 "0 1~' "1 = "n — 

z[l — z) z z — 1 l + oo 

A direct calculation shows that 

A[y]' (-1/ao) - ai(2 + a„)A[y] (-1/ao) =0, oi = -^ 

Hence 

A[f(^)]=£«(-ai(2 + ao);-l/ao) 
The latter is an Xi flag, and Xi flags are 1-step. Therefore, the desired conclusion follows by Lemma 
HH D 

Lemma 4.8. The flag £q^ is a 1-step flag. 

Proof. Define 



A[y] := ^ 



z(z-l) 
Using (|18cp . a direct calculation shows that 

where the last equality should be understood as an equality between polynomial flags. The desired 
conclusion follows by Lemma HH] D 

Lemma 4.9. The flag £{2 is a 2-step flag. 

Proof. The degree regular basis is shown in (|18dl) . In particular, note that 

C/i = span{2z — 1}. 

Define 

aiW[y,2z-l] 

My] — r. — ^ — 

z(l-z) 

A direct calculation shows that 

AM'(l/2)=0. 
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Hence 

The latter is an Xi flag, and Xi flags are 1-step. Therefore, the desired conclusion follows bv 14.41 D 

Lemma 4.10. The flag ^23 (a) is a 2-step flag. 
Proof. Define the operator 



A[y]:^{y'^ay)/z + Ky', K=Va^±3 

Applying A to the degree regular basis shown in (j28ep gives a flag with a stable degree sequence of 1, 2, . . .. 
Imposing 

y'(0) = 02/(0), y"'{0) = 3ay"{0) ± 6ay(0), 
a direct calculation shows that 

A[ynO) = {a + K)A[y][0]. 

Since the former conditions defines f (^^) and since the latter conditions defines 8^^'>{a + K;0) (see ^ for 
the definition), it follows that 

A[£(2'')]c£(i)(a + X;0) 

Next, define 

B[y] := z(l - Kz)y' - (3 + (a - 2K)z)y 
If we suppose that 

y'{0) = {a + K)y{0) 
then by direct calculation, 

B[y]'iO) = aB[ym, B[y]"'iO) = 3ay"i0) + 6ay(0) 

Therefore, 

B[£-(i)(a + if;0)]cf(2fc). 

Next, observe that A lowers degree by 1 and that B raises degree by 1. Hence BA and AB do not 
raise degree and they preserve their respective flags. Since £{a + K; 0) is a 1-step flag fTheorem l3.1l and 

f(36 
^23 



Lemma [4T]l it follows that £23 is a 2-step flag. D 



Lemma 4.11. The flag £23^' (a) is a 2-step flag. 

Proof. The argument is the same as for the proof of Lemma |4.10| but with the following operators: 

A[y]:=^— + ^y 

B[y] := z(l + (a - l)z)y' - (3 -f (2a - l)z)y. 

We then have 

A[£^^'^]c£^°\l;0), B[f(°)(l;0)] cf('"^ 

D 

Lemma 4.12. The flags ^13° , ^ 03° , ^12" 1 1^02" '^^^ '^^^ 2-step flags. 



Proof. By Proposition 13.41 all of the above flags are various specializations of 

£*^^^(aoi,0,a23;0) = spanjl + aoiz,z^ + a2s.z^ , z'^ , z^ , . . .}. 

Hence, it suffices to prove the assertion for this general case. Equivalently, the above flag consists of 
polynomials satisfying 

y'(0) = aoi2/(0), y"'[Q) - 3a232/"(0) (38) 

Consider the operator 

^[y\ := ^ aoi?/ 

and note that 

^[aoiZ + 1] =0. 
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Next, observe that 



A[y]'{z) - i(aoi + 3a23)A[2/](z) = (aoiy(O) - y'(0)) (^^ + 



(floi + 3a23)/2 



+ \{y"'{0)-ia^,y"m + O{z) 

Hence, if y{z) satisfies dSH]), tlien A[y\ e f (^^(aoi + 3a23)/2; 0). 
At this point, let us suppose that aoi 7^ and note that 

A\ VI \ A\ M \ i^ ^ \ /^"Ol \ + aQlZ , I , 

A[y\ [z) - aoiA[y\{z) = [I + aoiz) -^y ^ V + -V 

V z^ z^ z 



Hence A[y] e .^'^'^•'(aoi; --1/aoi) for ah polynomials y{z). Together, the above calculations demonstrate 
that if aoi 7^ 0, then 

^[^(2a)] ^ £(ii)((ao, + 3a23)/2, aoi; 0, -1/aoi). 

Hence, by Lemma H^ the flags ^13" , 1^12'' are Darboux connected to the flag above. We already showed 
that £(^^) is a f-step flag, so this concludes the proof for the case aoi 9^ 0. 
Finally, let us consider the case aoi =0. In this case, 

A[y] = ^, 

z 
£(2)(o, 0, 023; 0) = span{l, z^ + a23z3, z\ z\ . . .} 

A[£(2)(0, 0, a23; 0)] = span{2 + 8023^, ^^ 2^ • • •} 
= £:(!) (3a23/2;0) 

By Lemma l4Tn the latter is a I-step flag. Since A[l] = 0, applying Lemma \4A\ shows that Sq^"" ,£02^ are 
both 2-step flags. D 



Proof of TheoreTn \4.S\ There are two basic mechanisms which we use to give the proof of the conjecture 
for X2 and Xi operators. The first mechanism is that of dimensional exhaustion, and is utilized in Lemma 
14.21 and in Lemma 14.51 This mechanism is used to prove the conjecture for Xi flags (Lemma 14. 2p and 
also used in the proof of Lemmas 14. 6[ 14.101 and 14.111 All these cases require that we exhibit both an A 
operator, which relates the given flag U to a. "simpler" flag U, and a B operator that relates U back to U. 
The other basic argument is conceptually related to state-deleting transformations in quantum me- 
chanics. Here it suffices to show that a 1st order operator that annihilates Ui maps the given flag U to 
a simpler flag U and to have in hand a second order operator that preserves the given U. This is the 
argument of Lemma l4.4l This argument is utilized in Lemmas l4.71l4.8l and l4.9l4.12l Taken together, these 
Lemmas cover the cases of all possible Xi and X2 flags and the operators that preserve them. D 

5. POLYNOMIAL StURM-LIOUVILLE PROBLEMS AND DARBOUX TRANSFORMATIONS 

Our main goal is to complete the classification of X2 OPS and what remains to do is to select from 
all the X2 operators given in Section [3] for each X2 flag, those that give rise to a well defined Sturm 
Liouville problem. For this reason, in this Section we need to review some preliminary results from the 
theory of Sturm Liouville problems. We will also provide the main definitions and properties of algebraic 
Darboux transformations for second order differential operators. We emphasize that by construction these 
transformations will map an SL-OPS into an SL-OPS. 

5.1. Orthogonal polynomials on the real line defined by a Sturm-Liouville problem. Every 
second-order eigenvalue equation 

T[y] := p{z)y" + q(z)y' + r{z)y = \y 

can be put into formal Sturm-Liouville form 

-{Py'y + Ry ^ -\Wy 
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where 



P{z) — exp / q/pdz 

Wiz) = (P/p)(z), 
R{z) = -irW){z), 



(39) 

(40) 
(41) 



With the above definitions, the operator T[y] is formally self-adjoint with respect to the weight W{z)dz 
in the sense that Green's formula, below, holds: 



j T[f]g Wdz ~ J T[g]f Wdz = P{f'g - Jg') 



(42) 



If the operator T[y\ has infinitely many polynomial eigenf unctions, and if an interval of orthogonahty can 
be appropriately chosen so that W{z)dz has finite moments and the right-hand side of (l42t vanishes for 
polynomials f{z),g{z), then the eigenpolynomials of T[y] constitute an SL-OPS. 

By direct inspection, every X2 operator listed in Propositions 13.21 and 13.51 has the form 

T[y] := p{z){y" - 2(logO') + qiz)y' + r{z)y 

where p{z) is a quadratic polynomial, q{z) is a linear form, ^(z) is either z(z— 1) or z and r{z) is a rational 
function with ^(z) in the denominator. Applying an affine change of variable, 

z — ax + b 

the coefficients p{z) and q{z) can be put into a normal form. There are five classes of these normal forms, 
which we display in Table 1 together with the interval of orthogonality and the weight defined by (p9| 

-(SOD. 







Table 1. 






p{x) 


q{x) 


W{x) 


I 


OPS family 


1 


-2a; 


l(x)- 


(— CX3, 00) 


Hermite 


X 


a -1- 1 — x 


ax)- 


(0,00) 


Laguerre 


l-x^ 


P-a-{2 + a + P)x 


(i-xni+xf 

iixy 


(-1,1) 


Jacobi 



X 


2(a;±l) 


ax)- 


n/a 


Bessel 


\ + x^ 


a + 2{p + l)x 


ax)- 


n/a 


twisted Jacobi 



Just as in the analysis of classical orthogonal polynomial systems [23] the Bessel and twisted Jacobi 
cases can be excluded because it is not possible to choose a interval of orthogonality that satisfies the 
finite-moment condition. Therefore the search for X2 orthogonal polynomial systems narrows to the 
first 3 cases. In each case, the requirement is that S,{z) have no zeros on the corresponding interval of 
orthogonality. For the Laguerre subcase, there is the additional constraint that a > — 1. For the Jacobi 
subcase, the constraint is that a,/3 > —1. 



5.2. Factorization and orthogonal polynomials. Consider two differential operators: 

T[y] ^py" + qy' + ry, 
T[y] = py" + qy' + fy, 
related by a factorization (pO)) (ISTI) . Let us write 

A[y] = Ky' - wy), 

B[y] == b{y' -wy), 
where p{z),q{z),r{z),b{z),w{z),b{z),w{z) are all rational functions. We will refer to 

(f>{z) = exp / wdz, w = 4)' / (f) 



(43) 
(44) 

(45) 
(46) 

(47) 
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as a quasi-rational factorization eigenfunction and to b(z) as the factorization gauge. The reason for the 
above terminology is as follows. By (PO]) . 

T[4>] = Ao0; (48) 

hence the term factorization eigenfunction. Next, consider two factorization gauges bi{z),b2{z) and let 
Ti[y],T2[y] be the corresponding partner operators. Then, 

T2 = fi~^TiiJ,, where ^(z) = bi{z)/b2{z). 

Therefore, the choice of b{z) determines the gauge of the partner operator. This is why we refer to b{z) 
as the factorization gauge. 

Proposition 5.1. Let T[y] be a second order rational operator that preserves a polynomial flag. Let 
4>{z) be a quasi-rational factorization eigenfunction with eigenvalue Xq. Then, there exists a rational 
factorization (1301) such that the partner operator preserves a primitive polynomial flag. 

Proof. Let w{z) = (j)' (z) / (j){z) and let b[z) be an as yet unspecified rational function. Set 

w^~w-q/p + b'/b, (49) 

b^p/b, (50) 

and let A[y],i3[2/] be as shown in in (j45|) (|46|). An elementary calculation shows that (|30|) holds. Let 
yi,2/2,--- be a degree- regular basis of the flag preserved by T. We require that the flag spanned by 
A[yj] be polynomial and primitive (no common factors). Observe that if we take b{z) to be the reduced 
denominator oi w{z), then A[yj] is a polynomial for all j. However, this does not guarantee that A[yj] is 
free of a common factor. That is indeed a stronger condition which in fact flxes the gauge b{z) up to a 
choice of scalar multiple. Finally, the intertwining relation 

TA = AT (51) 

implies that A[yj] are eigenpolynomials of the partner T. D 

In the preceding subsection, we showed that a second-order operator T[y] is formally self-adjoint relative 
to a weight W defined by ([M]) (HH]). The following Proposition describes the effect of a factorization 
transformation on the corresponding factorization function and the weight. 

Proposition 5.2. Suppose that rational operators 

T[y] = py" + qy' + ry, f[y]^ py" + qy' + ry 

are related by a rational factorization with factorization eigenfunction (t>{z) and factorization gauge b{z), 
Then the dual factorization gauge, factorization eigenfunction and weight function are given by 

bb^p (52) 

W/b = W/b, (53) 

1$ = 1/{W(I3) (54) 

Proof. Equation (|52p follows immediately from (l45t (|46l) ((30)) . From there, equation ([3T|) implies that 

w -f w = -q/p + b'/b = -q/p + b'/b. (55) 

Hence, 

q = q+p' ~2pb'/b. (56) 

From here, ^^ follows by equations ^^ (gHl). Equation dM]) follows from (gT]). D 

The dual weights W, W allow us to interpret the intertwining operators A[y], B[y] in terms of a formally 
adjoint relation 

J A[f]g Wdx + J B[g]f Wdx - iP/b)fg (57) 

If the right hand side vanishes on an appropriately chosen interval of orthogonality, and if the partner 
operators T, T both admit an infinite sequence of eigenpolynomials, then the operators T and T and their 
corresponding eigenfunctions are related by a 1-step Darboux transformation. 
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The dual factorization functions (j), </> allow us to express the adjoint intertwiners as Wronskians: 

A[y]^b(t>-'W[4>,y] 

B[y] = br'W[4',y]. 

In Theoreni l4.2l of Section |4l we established that every X2-operator is Darboux-connected to a classical 
operator and that the requisite intertwiners also connect the corresponding exceptional flag with the 
standard polynomial flag. Theorem 11.11 follows as an immediate corollary. In light of the above remarks, 
it is convenient to give the connecting intertwiners as Wronskians of factorizing functions of the classical 
operators. Therefore, before turning to the exhaustive classification, we must review the possible quasi- 
rational factorizing functions for the classical operators. 

5.3. The X2 Hermite polynomials. The classical Hermite orthogonal polynomials are orthogonal rel- 
ative to the weight 

The nth Hermite polynomial H„{x) satisfies the differential equation 

n[H„] = -2nH„ 

where 

H[y]=y"-2xy' 

The exhaustive classification of the X2 polynomials confirms the factorization conjecture. This means 
that all X2 Hermite polynomials are given as Wronskians of the classical polynomials together with fixed 
quasi-rational factorization eigenfunction of the classical Hermite operator 'H[y]. These quasi-rational 
eigenfunctions are listed below: 

^i'\x) = Hr^ix), H[7A(i)] = -2nV'(i) (58) 

We will use Hn{x) to denote the X2 Hermite polynomials, where the degree index n skips exactly two 
values. These exceptional Hermite polynomials are orthogonal relative to a weight of the form 

W{x]a,fi) 



a^) 



where the denominator ^(x) is a quadratic polynomial. Consequently, the Hn{x) are eigenpolynomials of 
an operator of the form 

n[y]:='H[y]-2{\ogO'y' + r{x)y, 

where r(x) is rational in x and where the prime denotes a derivative with respect to x. In order for the 
weight to be non-singualr, the quadratic ^(x) must have imaginary roots. Also, as we show below, the 
rational term r{x) = always vanishes. This is established on a case-by-case basis, and has no apriori 
explanation. 

5.4. X2-Laguerre polynomials. The classical Laguerre weight is 

Wc,{x) ^ e-^x" 
The classical Laguerre operator is 

Ca[y\ ■■= xy" + (a + 1 - x)y' 
The quasi-rational eigenfunctions of this operator are 

</)li) (x; a) - i(f ) (x) £„ [</>,W] = -n</>,W (60) 

</.(2)(x;a) =.T-"4-")(x) £„[</)(2)] = (a-n)0(2) (61) 

0(f)(a;;a) = e-4")(-x) CM'^] = ia + n + l)<lyl^^ (62) 

^';^\x;a) = e^x~'^Ll-'^\-x) C^^] = (n + 1)^('^) (63) 
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In confirmation of the factorization conjecture, all X2 Laguerre polynomials are given as first and 
second-order Wronskians of the classical Laguerres and the above factorization functions. The X2 poly- 
nomials themselves will be denoted by Ln the range of n omits exactly two degrees. In all cases, the 
Ln are orthogonal relative to a weight of the form 

mx;a):=|^, (64) 

where the denominator ^(x; a) is a quadratic polynomial in x. The parameter a has to be restricted so 
that ^(x;a) has no zeros in the interval of orthogonality x G (0,oo). The exceptional polynomials I/„ 
arise as eigenpolynomials of a second order operator 

£„ [y] = xy" + (1 + a - x)y' - 2{\og^)'y' + r{x] a)y 

where r(x] a) is a rational function in x which will be adjusted so that, in all cases, 

£„[l1")] = -nil") 

5.5. The X2 Jacobi polynomials. The classical Jacobi OPs are orthogonal relative to the weight 

W{x;a,P)^{l-x)°'{l + xf, a,/3>-l. 
The nth Jacobi polynomial Pn ' [x) satisfies the differential equation 

ra,/5[Pi"'^'] = -n{n + a + P + l)Pi"'« 

where 

Tc,ji[y] = (1 - x')y" + {l3-a-{2 + a + P)x)y' 

The exhaustive classification of the X2 polynomials confirms the factorization conjecture. This means 
that all X2 Jacobi polynomials are given as Wronskians of the classical polynomials together with fixed 
quasi-rational factorization eigenfunction of the classical Jacobi operator 7^,/3[2/]. These quasi-rational 
eigcnfunctions are listed below: 

(/.ii'(a;;a,/3)=p("'^'(a;), Ti^'^^^] ^ -n{n + a + p + 1)^'^^^ (65) 

(/.i2)(a;;a,/3) = (l+a;)-^p("^-^)(x), r[(/)(2'] - (/3 - n)(n + a + l)(/)(2) (66) 

^'i\x;a,l3) = il-x)-"P^-"-^^Hx), r[(/)(^'] = (a - n)(n + ;3 + l)(/.(3) (67) 

4'*) (a;; a, /?) = (!- x)-"(l + a;)-^p(-"'-«(a;), r[(/)(^'] = (n + l)(a + /? - n)^/.^^) (68) 

We will use P„ (x) to denote the X2 Jacobi polynomials, where the degree index n skips exactly two 
values. These exceptional Jacobi polynomials are orthogonal relative to a weight of the form 

W{x;a,l3) = — ^ 

where the denominator ^(x;a,/3) is a quadratic polynomial. Consequently, the Pn {x) are eigenpoly- 
nomials of an operator of the form 

ra,,3M :=ra,^[y] -2(l-a;2)(loge)V + 7-(a;;a,/3)2/, 

where r(x; a, /3) is rational in x and where the prime denotes a derivative with respect to x. The parameters 
a,j3 > —1 are so restricted in order to have finite moments of all orders. Additional restrictions must be 
imposed on a, j3 to ensure that ^(a;; a, j3) has no zeros in the interval of orthogonality x G (—1, 1). 

6. Classification of codimension 2 XOPs 

The main result of this section is a complete list of X2 orthogonal polynomial systems together with 
the intertwining operators that connect them to the classical families of Hermite, Laguerre and Jacobi. 
The classification is summarized in the following. 

Theorem 6.1. Up to a real affine transformation of the independent variable, all X2 orthogonal polyno- 
mial systems are gathered in the following table: 
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Table 2. Classification of X2 orthogonal polynomial systems 





c(ll) 
'-'23 


c(ll) 
<^13 


c(ll) 
"--03 


c(2a) 
'-■13 


M2a) 
"--03 


Hermite 






1-step 
^6.1.^ 






Laguerre 


1-step 
^6.2.1\ 


2-step 
i6'.i^.i>| 


1-step 
^6.2.3\ 


2-step 
i\6.2.5\ 


2-step 
^6.2.6\ 


Jacobi 


1-step 
^6.3.1\ 


2-step 


1-step 
fJ6'.5..V| 


2-step 
i6'.5.5| 


2-step 



In Table 2 we find the classification of X2 orthogonal polynomial systems. In each cell we give the 
number of iterated Darboux transformations to obtain these families from a classical OPS, and we specify 
the subsection where each family is described. Empty cells mean that an OPS of that type does not exist 
for the given flag, and the same is true for all the other X2 flags not included in the table. The cells 
marked in bold correspond to X2-OPS previously known in the literature, while all other cases are new. 

In the rest of this section we will select the X2 operators for each of the X2 flags in Section [3] that 
can be transformed into a well defined Sturm Liouville problem of Hermite, Laguerre or Jacobi type. We 
allow affine changes of variables and basically we need to transform the leading order of the X2 operator 
into 1, a; or 1 — a;^ and verify that the weight is non-singular in the corresponding interval and it has 
well defined moments of all orders. This will exclude many cases and it will impose constraints on the 
remaining free parameters for the cases that survive. 

6.1. X2-Hermite OPS. 

6.1.1. No Hermite polynomials for the 2-pole flag ^23 • The leading order coefficient in (|19a|) . is 

1 a^ 

--z^(ao - ai)(ao - ai + 4) - z{aoai - oq - of + 3ai) - ^ + «i 

We require the coefficient of z^ to vanish. Setting ai — qq transforms the above into 

-ao(2z + ao/2- 1) 

Setting ai = ao + 4 gives 

(ao + 2)(2z-ao/2 + 2) 
In other case, it is impossible to obtain a Hermite-like operator. 



6.1.2. No Hermite polynomials for the 2-pole flag £\^ . The leading order coefficient in (J19bp . is 

-,2 / i- 



-(co + ci)y+co ( z ^ 

It is not possible to specialize co, ci so that the above polynomial reduces to a constant. 

6.1.3. 1-step Hermite polynomials that span the 2-pole flag E^^ . Setting ao,Q:i — —1/2, go = 1 hi (|19cl) 
and applying the change of varibles 

z = ^/^/2x + 1/2 
gives a Hermite-type operator 

U[y]:=y" -2xy-2{\ogO'y' 
where 

i{x) = l + 2x' = -\H2{ix). 
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The adjoint intertwiners and the exceptional polynomials are shown below: 

B[y]=e--V[V'(2),2/] (69) 

m = {-^ (70) 

-ffo = 1 (71) 

i?„=B[i/„_3], n = 3,4,5,... (72) 

HlHr,] = -2nHn (73) 

A[7J„] == 4ni/„^3, n = 0,3,4,5,... (74) 

The above polynomials are related to the CPRS exactly-solvable potential [31 [7] and constitute the 
codimension-2 instance of the modified Hermite polynomials introduced in [3]. This family was also 
described independently in |S^ for arbitrary codimension. 

6.1.4. No Hermite polynomials for the 2-pole flag £^2 ■ By inspection of (J19dp . a Hermite-type operator 
requires 

ao = ai = -go 7^ 

Applying a change of variable 

z = ax 
yields the weight 

W^(a:) = e-2-'-'(l-4aV)2 

To have a real weight requires a to be either real, or purely imaginary. In the first, case, the weight is 
singular; in the latter case there are no singularities but the finite-moment condition is violated. 

6.1.5. No Hermite-type polynomials for the 1-pole flags £(2a)^^(26) ^^^ E^^^i. A real valued operator and 
weight requires the unique pole to be real. However, a Hermite-type weight requires the entire real line 
as the interval of orthogonality. Therefore, even if Hermite type weights of the form 

do exist, since h is real, the resulting weight is singular. 
6.2. Xa-Laguerre OPS. 

6.2.1. 1-step Laguerre polynomials that span the 2-pole flag ^23 • By direct inspection of (|19a|) a Laguerre- 
type operator requires either ai = ao -I- 4 or ai — a^. We consider these two cases in turn 
(I) Imposing ai = oq -l- 4 in (|19ap . making an affine change of variable 

x = (ao-l-2)(2-ao/4-l), 

and setting 

a = ao(4-t-ao)/4 
gives the operator 

CM := xy" + (1 + a - x)y' - 2{\ogC)' {xy' + ay) 

where 

^{x; a) = 4"~^'(-a;) = {x^ + 2{a + l)x + a[a + l))/2 
and the prime symbol denotes the derivative with respect to x. We impose a > in order to 
avoid positive zeros of ^(x; a). The resulting orthogonal polynomials are codimension-2 instances 
of the type I exceptional Laguerre polynomials il4 t i28 ) . The corresponding polynomials and the 
adjoint intertwining relation are shown below: 

A[y]:^x^+^W[x~^,y]/i{x-a) (75) 

B[y]:^e--W[4>f\x;a-l),y] (76) 

L'CHx) = B[L^:i,\ n = 2,3,4,... (77) 

4L(")] = (a + n)Lt-'^. (78) 
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(II) Imposing ai ~ oq in (J19ap . making an affine change of variable 

X = ao(4z-2 + ao)/4, 

and setting 

a = all A - 1 

gives the operator 

ta[y\ := xy" + (1 + a - x)y' - 2x (log^)' {y' - y) 

where 

i{x] a) = 4~"~^^(a;) = [x"^ + 2(a - l)x + o? - a)/2, a > 1 

The resulting orthogonal polynomials are codimension-2 instances of the type II exceptional La- 
guerre polynomials [T4lf28] . The definition of these polynomials and the adjoint differential relation 
are shown below 

A[y]:^^^W[e^,y] (79) 

B[y]:=x''+^W[<j>^i\x;a + l),y] (80) 

U^^ = B[lI^+,% n = 2,3,4,5,... (81) 

A[L(r)] = (3-a-n)LtY^ (82) 

6.2.2. 2-step Laguerre polynomials that span the 2-poles flag £[^ . By direct inspection of (J19bp . a 
Laguerre-type operator requires cq — l,ci — 0. Applying the affine transformation 

^=(z-l/2) °°^^ + "°^ 
ao + 1 

and setting 

_ gg + 2ao + 2 
2(ao + l) 
gives the operator 

Co.[y] := xy" + (1 + a - x)y' - 2a:(logO'y' + ^(" " ^H" + ^ - ^) y 

$(x;a) 

where 

^(x; a) = .T^ + 1 - a^ = e-^^x^+^Wi^^^^ (2:; a), 0^ (x; a)], |a| < 1 

The adjoint intertwiners and the exceptional polynomials are: 

B^ [y] := ie-2.^2+oy^;[^(3) (^. ^^^ ^(4) ^^, ^^^ ^j ^33 

a 
l["^ ■.^L["\~x)^x + a + l (84 

L(r):=i?„[4"J3], n = 3,4,5,... (85 

My] -^7 ^W[x~"(a:;-a + l),.T + a + l,2/] (86 

A„[L(,")]=-(n-l)(a + n-l)4"J3, n = 1, 3, 4, 5, . . . (87 

Note: for a = 0, the above definitions have to be treated as a limit process. A straightforward calculation 
shows that 

BoM = -x{\ + x^)y" + (2x^ +x^ + 2x~ \)y' - {x^ +x^ + 2x- 2)y 
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6.2.3. 1-step Laguerre polynomials that span the 2-poles flag £q^ . Inspection of (|19cp reveals that a 
Laguerre-type operator requires 

qo + co + ci = 1 
Since we are free to scale the operator, no generality is lost by imposing cq — ci = 1, which gives us 

co = (l-9o)/2, ci = -(l + go)/2 

Applying the afhne change of variables 

X = go(l -9o - 2z) 

and setting 

a = 1 - 90 
gives the operator 

4 [y] := xy" + (1 + a - x)y' - 2a;(log ^'y' (89) 

where 

^{x; a) = 4"""^ V-^) = (x^ + 2(1 - a)x + a^ - a)/2, (90) 

and where 

a e (-1,0)U (l,oo) 
in order to avoid positive zeros in ^{x; a) and to have finite moments. The corresponding exceptional 
polynomials and intertwiners are shown below: 

B[y] :=e-"a;2+"W[4^^(:E;l+a),2/] (91) 

A[y] = 77^ (92) 

l[-\x) = 1 (93) 

L(^\x) ^ B[l},:+^\ n = 3,4,5,... (94) 

A[4")] = nLi"+/\ n = 0,3,4,5,... (95) 

6.2.4. No Laguerre polynomials for the 2-poles flag £\^2 ■ By inspection of (J19dp . go = cq + ci. Without 
loss of generality, 

Co — ci = 1, Co + ci = a 
where a is a new operator parameter. Making the affine change of variables 

X = a((l + a) - 2z) 
gives the weight 

^"^^^ " {x-a^-aY{x~a-^+aY 
In order to have a real weight we need a to be either real or pure imaginary. In the first case, the 
denominator will have a positive zero; the weight is singular. In the former case, the finite moment 
condition is violated. Therefore, there are no X2 polynomials that span this flag. 

6.2.5. 1-step Laguerre polynomials for the 1-pole flag £1^ ■ We refer to the S^"^' flags and the correspond- 
ing OPS as 1-pole because the weight function has one pole, unlike the 2-poles present in the weight 
functions of the £^(^^-' families. This pole in the weight has higher multiplicity. 

By direct inspection of (I29ap . a Laguerre-type operator requires either a = 1/3, or a = 3. Setting 
a = 1/3, making the change of variables x = z + 3/4 yields a singular weight, namely 

g-x^-l/4 



{Ax - 3)4 

Setting = 3 and making the change of variables 

X = 3z - 3/4 

gives the operator 

£M:=.Ty" + (5/4-x)y'-^^^ 

X -\- 3/4 
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and the weight 

p-x IjA 

which is both non-singular and has finite moments of all orders. The remarkable feature of this weight is 
that it has a fourth order pole, unlike the two second order poles of the previously discussed X2 families. 
The adjoint intertwiners and the exceptional polynomials for this weight are shown below: 

— 1x 9/4 

B\V\ := ■^^|^W[<?ii')(a:;l/4),4^)(x;l/4),y] (96) 

Li(j:) :=x + 15/4 (97) 

K(x) -^ B{L^IL% n = 3,4,5,... (98) 

a;9/4 

(x + 3/4) 
^[in] = ^(n-l)(4n+l)L,t/^), n = 1,3,4,5,... (100) 

6.2.6. 2-step Laguerre polynomials for the 1-pole flag £q^ ■ By inspection of (j29b|) . a Laguerre-type op- 
erator requires go = 3. Making the affine change of variable 

a:=-(2z-l) 
4 

gives the operator 

t[y] := xy" + (3/4 - x) y' - _^ ^ ^ 

and the weight 

„-x -1/4 

The adjoint intertwiners and the exceptional polynomials are shown below: 

— 2x 7/4 

^[y] ■= ^-^W[#(^;-l/4),0f^(x;-l/4),2/] (101) 

X -I- 3/4 

Lo = 1 (102) 

U, := B[L\;}i% n = 3,4,5,... (103) 

7/4 

^[y] '-^ r^^TTT^s W[l' ^'^'(^ + 15/4), y] (104) 

(x -I- 3/4) 

A[Lr.] = ^n(5 ~ 4n)4-V'^ n = 0, 3,4, 5, . . . (105) 

6.2.7. No Laguerre polynomials for the 1-pole flags 1^02" , ^12" , '^23 ^^^ ^2'3^ ■ Setting po = and apply- 
ing an affine transformation, the operator (j29cp yields a singular Laguerre-type weight 



W{x) 



g-:r 3,1/4 



[Ax - 3)4 
By direct inspection of (I29dp . (|29ep (|29f[) . the operators in question do not admit a Laguerre form. 

6.3. X2 -Jacobi OPS. 

6.3.1. 1-step Jacobi polynomials that span the 2-pole flag £23 • The quadratic coefficient of y" in (|19al) 
factors as 

- -(ai - ao)(ai - ao - 4)(z - zi)(z - Z2) (106) 

where 

"1 ai-2 

2^1 = 7, ^2 = T (107) 

tti ^ flo ^ 4 oi — ao — 4 
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We seek an affinc change of variable that transforms this quadratic into l~x^. There are two possibihties 
according to which root is sent to +1 or —1. However, since the two resulting families are related by an 
af&ne change of variable, it suffices to consider just one such transformation. Employing the transformation 

z = |(x + l)-|(x-l) 



setting 



2(Z1-1)Z1(2Z2-1) _ 2(2zi - 1)Z2(Z2 - 1) 

Zi- Z2 ' Zl - Z2 

.(ll)r 



and adding a constant term, transforms T2^ [y] into the operator 

faAy] = To^Ay] - 2(logO'((l - x^)y' + /?(! - x)y) + 2(a ~ /3 - l)y 
where 

ax;a,P)^Pt''-'^^-'\x) (108) 

= J(^"2^^)(--l)^ + ^(/3-« + l)(l-«)(x-l)+Q (109) 

In this way, we have arrived at the codimension-2 instance of the exceptional Jacobi-type polynomials 
introduced by Odake and Sasaki [I6 1 I28 ] . 

We require that ^{x;a,l3) have no zeros in the interval of orthogonality x € [—1, 1]. The above affine 
transformation maps —1, 1 to the roots of ^{x) and maps 



Zl 



- ± -^/-a{l + a + b)/b, a^a-1, b=-/3-l (110) 



2 2 



^2 = ^T^V-&(l + a + &)/a (111) 

to ±1. Therefore, an equivalent condition is that zi,Z2 are either complex- valued or lie in the interval 
(0, 1). The solutions to this constraint in the (a, 6) plane are the disjoint union of the following regions: 
(i) a,b > 0; (ii) a > 0, 6 < — 1; (iii) a < —1, 5 > 0; (iv) — 1 < a, 6 < 0. Finite moments require a,/3 > —1. 
Therefore, in the final analysis, we have two classes orthogonal polynomials with a non-singular weight 
and finite moments: a > —1,/? > and < a < 1, — 1 < /3 < 0; c.f.. Proposition 4.5 of [16] . 
The exceptional polynomials and the adjoint intertwiners are shown below: 

A[y]:=^-^±^^W[{l + x)-^,y] (112) 

B[y]:={l-xr+^W[4^^^\x;a + l,^~l),y] (113) 

Pi"-« = B[Ptt'-^^''^] (114) 

fa,0=BA+{2 + l3){a-l) (115) 

Tc,+i^l3^i = AB + {2 + I3){a - I) (116) 

r[P„] = -(n-2)(n-l + a + /?)A. (117) 

^[P("-'^)] = -(a + n- 3)(/3 + n)Pi'^_1,'-'-'^ (118) 

6.3.2. 2-step Jacobi polynomials that span the 2-pole flag £i^ . The quadratic coefficient of y" in (|19bp 
factors as 



-{{R + l)z - l)((i? - l)z + 1), where R = J -- 
2 V Co 



Employing the affine transformation 
and setting 



_ i?a;-|-l 
1 flo R 



1-R 1-R {l + ao)il-R) 



(119) 



^ " TTi? ^ TTi? ^ {l + ao){l + R) ^^^^^ 
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transforms the operator T^g into 

^ ri ^ ri „. 2^n .n/, 8(a - 1) (/3 - l)P|"^'^Va;) 

r„,M2/] = ^a,M2/] - 2(1 - x2)(iogoy - ^ — t(i^^if — y 

where 

C(a;;a,/3) = {x^ + l){a^ - p'^) + 2x{a^ + p'^ - 2) (121) 

4ao(2 + ao)(l + gp - R){1 + ap + R) ^ 
= (l + ao)^(fi^-l)^ ^^ + "^^^"^^ + '^ (^''^ 

For a real, non-singular weight, we require i? = e'*, t e K to be a unit-length complex number. A direct 
calculation shows that 



a2 _ ^2 ^2 _ ^2 



1 _ a^+p^-2 2V(a2 - l)(/32 _ j) 
i? ^ a2 - ;32 T a2 - /32 

Therefore, the parameters a, /3 must satisfy 

-1 < a < 1, (3 > 1, or a > 1, -1 < ^ < 1 

The corresponding exceptional polynomial and the adoint intertwiners are shown below 

B[y]:=^ p W[cl)^^\x;a + 2,f3),cf,['\x;a + 2,l3),y] (124) 

p("'« = p|-"-2.^) (125) 



p(a:/3) ^ 5 



p("+2,/ 



71 = 3,4,5,... (126) 



f[Pn] = -n(n - 3 + a + /?)P„ (127) 

^[^i"'''^] = ^("-l)(^ + «-2)(n + /3-l)(7i + a + /?-2)Pi^f'''\ n = 1, 3,4, 5, . . . (128) 

As above, for the case oi P — 0, the definitions above must be treated as a limit. 

6.3.3. 1-step Jacobi polynomials that span the 2-pole flag £q^ . The quadratic coefficient of y" in (J19cl) 
factors as 

- (go + CO + ci) y + ^ + CO (^z - = -^^{z - ^i)(^ - ^2) (129) 

where, 

co = l, ci= fl--Vl--V qo = -2 + - + - (130) 

\ ZlJ \ Z2J Zl Z2 

Note that no generality is lost by scaling cq = 1 because, if Cq = 0, then the operator does not have a pole 
at z = 0. Employing the affine transformation 

zi{x + 1) - 22(1 — x) 

' = 2 

and setting 

2(zi - 1)Z1(2Z2 - 1) . 2(2zi - 1)Z2(Z2 - 1) 

a = , p - 



Zl - Z2 Zl - Z2 

and adding a constant term, transforms 2ziZ2T'23 ii^to the operator 

r„,^[y] = r„,/3[?/]-2(iogO'(i-x2)y' 
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where 

ax-.a,P)=Pt''-^-^-^\x) (131) 

1 /2 - /3 - a^ 



(x - If + 2(1 - /? - «)(1 - «)(^ - 1) + (^2 j (132) 

We require that ^(x; a, /?) have no zeros in the interval of orthogonality x E [—1, 1]. The above afHne 
transformation maps —1,1 to the roots of ^(x) and maps 

zi = -± 2 V-a(l + a + b)/b, a = a- 1,6 = ^-1 (133) 

22 = ^ T iV-fe(l + a + &)/a (134) 

to ±1. Therefore, an equivalent condition is that zi, Z2 are either complex-valued or lie in the interval (0, 1). 
This constraint, toghether with the finite moments constraint, gives us 4 disjoint classes of acceptable 
parameter values: 

(i) a,/3> 1; 

(ii) 1< a< 3, -1 < ^ < 0, a + /3 < 2; 

(iii) 1< /3 < 3, -1< a < 0, a + /3 < 2; 

(iv) < a, /3 < 1 

The exceptional polynomials and the adjoint intertwiners are shown below: 

v' 
A\y] := —. T^-^-n (135) 

B[y] := (1 - x)2+"(l + xy+''W[4*\x; a + 1,(3 + l),y] (136) 

Po^"'^^ = 1 (137) 

Pi'^'^^ ^ B[P,tf''+% n^ 3,4,5,... (138) 

f[Pn]^-in-2){n-l + a + P)Pn (139) 

A[P^"'^)] = -n(a + n - 3)P^"+''''+'\ n = 0, 3, 4, 5, . . . (140) 

6.3.4. No Jacobi polynomials for the 2-pole flag £[2 ■ The quadratic coefficient of y" in (|19dl) factors as 
(CO + CI - ,0) ^ + (f - c,) - I = -^(. - .0(^ - ^.) (141) 

where, 

co-1, ci=fl--Vl--V go = 2---- + — (142) 

V ^1/ V ^2/ Zi Z2 Z1Z2 

Note that no generality is lost by scaling cq = 1 because, if cq — 0, then the operator does not have a pole 
at z = 0. Employing the affine transformation 

_ zi{x + 1) ~ 22(1 ~ x) 
^~ 2 

and setting 

2(zi-l)zi(2z2-l) ^ 2(2Z1-1)Z2(22-1) 
Zi - Z2 Zl - Z2 

gives a weight of the form 

Since 

zi = -±-\/a(l + a + fc)/^ a = a + 1,6 = ^ + 1 (143) 

1 1 



Z2 = - T -VKl + a + fe)/« (144) 

and since a, (3 > —1 is required for finite moments, the roots Zi, Z2 are real, and one of them lies outside 
the interval (0, 1). Therefore, if a,/3 > —1, the above weight must be singular on x £ (—1, 1). 
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6.3.5. 2-step Jacobi polynomials that span the 1-pole flag £i^ ■ The quadratic coefficient of y" in (|29ap 
factors as 

(1 - 3a) (3 - a) — + 2 (1 - a) z + l) = -((a - 3)z - 2)((3a - l)z - 2) 

In order to have a Jacobi-type operator, we require a 7^ 3, 1/3, —1; in the latter case we obtain a perfect 
square. Applying the affine transformation 

(.T + 1) X-l 



yields the operator 

where 
and where 



a — 3 3a — 1 

ra[y]:=r„,My]-4(i-2:2)(iogOV 



£,{x]a)' 
£,{x; a) = (1 + a)x + 2(a - 1) 



' fi 



a — 3 ' 3a — 1 

Just as for the Laguerre-type polynomials, the corresponding weight involves a 4th order pole: 

^{x;a)-^ 

In order to obtain a non-singular weight we must have a > 3 or a < 1/3. However, in order to have 
a, /3 > —1 (finite moments), we must restrict the latter condition toa<— 1/3, a^^— 1. The corresponding 
values of a,/? range from a > 2, 0</3<2in the former case, and 1/5 < a < 2, — 1 < /3 < 0, (a,/3) ^ 
(1/2, —1/2) in the latter case. Of course a, (3 are not independent, but rather are linked by the relation 

4a/3 + l3-a + 2 = 

The adjoint intertwiners and the exceptional polynomials for this flag and weight are shown below: 

B[v\--=-, L,KX . >V[c/)")(a:;a + 2,/?),4^^(x;a + 2,/3),z/] (145) 

a[a — 1)(1 + 6a)E,{x\ a) 

A[y]:^^^-^^^^f%W[{l + x)-PMl + a){x-l) + {a-l){2.a-l),y] (146) 

36(1 + 6a)^[x\ay 

Pi{x; a) = 2(1 + a){x - 1) + (a - l)(3a - 1) (147) 

P„ix;a):=B[P^lf], n- 3,4,5,..., (148) 

A[P„l]^{n-l){n~3 + a){n + f3){n-2 + a + l3)Pi"_:p, n = 1,3,4, 5, . . . (149) 

6.3.6. 2-step Jacobi polynomials that span the 1-pole flag ^03° ■ "^^^ quadratic coefficient of y" in (|29bl) 
factors as 

(3-,o):^-2z + l>-(^-^^)(^-^^) 
where 



-4±2yr+^ zi 

Zl,Z2 = , Z2 - 



90-3 22:1-1 

Applying the affine transformation 

zi{x + 1) (a; - l)z2 , 

rM:=ra,0[y]-4(l-x2)(loge)y 
^{x;zi) = (zi - l)x-\- zi 



yields the operator 

where 

and where 



3 3 

a = -zi - 1, /3=-Z2-l, 4Q:/3 + a + /3-2 = 
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Just as for the Laguerre-type polynomials, the corresponding weight involves a 4th order pole: 

In order to obtain a non-singular weight we require zi ^ Z2 to have the same sign. This implies that 
zi > 1/2, zi 7^ 1, which in turn implies that a, /3 > —1/4, a, /3 ^ 1/2 but are subject to the relation 

4a/3 + a + /3-2 = 

The finite moment condition is therefore automatically satisfied. The adjoint intertwiners and the excep- 
tional polynomials for this fiag and weight are shown below: 

(-1 _ ^^2Q+6('-| I ™^/3+2 
^M ■■- MP), ; W[4'^(x;a + 2,/3),0f)(a;;a + 2,/3),2/] (150) 

A[y] := (^_i)2^(^_2)2 >^[1' (1 + ^)"'^(1 + « + (^ - 1)/?(1 - 2«), 2/] (151) 

Po(a:^;^i) = l (152) 

h{x;zi):=B[Pi''_f'% n = 3,4,5,..., (153) 

A[F„] =n(n-2 + a)(n-l + ;3)(n-3-)-a-|-/3)Piil3"''^\ n = 0,3,4, 5, . . . (154) 

6.3.7. No Jacobi polynomials for the 1-pole flags ^02" j '^12° ' '^23 ''^23'^ • Setting 



-i±yr^^ 

2:1, 22 = 

Po 
and applying the afiine change of variables 

Zi{x +1) {x — 1)^2 
Z = . 

2 2 

transforms the operator in (|29cp into Jacobi form. The corresponding weight is 

(l-a;)°(l + a;)'^ 



zi 7^ Z2 



VF(a;;zi 
where 



(a;(zi + 1) + zi)4 



3 3 

a = -l + -zi, ^ = -l + -Z2 

A non-singular weight requires that zi , Z2 be real and have the same sign. Since 

-zi 

Zl = 

2zi + 1 

the only possibility is that zi,Z2 < —1/2. However, this means that a,/3 < —1, which violates the finite 
moments condition. 

By direct inspection of (|29dp (|29ep . a Jacobi- type operator must have a singularity at x = 0. The 
coefficient of y" in (j29fj) is a perfect square, which does not permit a Jacobi-type operator. 

7. Summary and outlook 

In the present paper we have given a classification of exceptional orthogonal polynomial systems of 
codimension two (X2-OPS). The classification includes all the cases previously known in codimension 
two plus some new examples of exceptional polynomials. Among the new families, the one-pole fiags are 
clearly special. Generically, the weight of a Xm-OPS is a rational modification of a classical weight with 
771 double poles, and this is the case for all the families known to date. The Jacobi and Laguerre OPS that 
span the E^ "> fiag have codimension two but only one pole in their weight, with quadruple multiplicity. 
They also have one less free parameter than the usual Laguerre and Jacobi families, i.e. no free parameters 
for the S^^"^ -Laguerre and just one free parameter for the f*^^"-'- Jacobi. The explanation for the presence 
of these exotic families is that generically they would belong to a higher-codimensional family, but that a 
careful tuning of the parameters can make the codimension drop by one and have two of the poles of the 
weight coalesce. Thus, the generic weight of an X^-OPS is a classical weight divided by the square of a 
certain degree m polynomial ^{x) with simple roots that lie outside the interval of orthogonality, but we 
know that degenerate cases are also possible. 
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We have also shown that every X2-OPS can be obtained from a classical OPS by a sequence of at 
most two Darboux transformations, and we conjecture this result to be true mutatis mutandis for any 
codimension m. Even if the conjecture could be proved to be true, the scheme of multiple step Darboux 
transformations is still very rich: there are four quasi-rational factorizing functions for the Laguerre and 
Jacobi families and two for the Hermite. The SL-OPS obtained by 1-step Darboux transformations have 
been studied in all cases, but multi-step Darboux transformations might mix factorizing functions of 
different kinds and all the possibilities have not yet been explored. It could also happen that even if the 
intermediate weights in a multi-step Darboux transformation are singular, the final weight will be regular. 
All cases when this happens have been studied for multi-step state-deleting Darboux transfomations 
in a more general Sturm-Liouville context (not necessarily polynomial) by Krein and Adler [T1I23). A 
generalization of Krein- Adlcr's Theorem to multi-step isospectral transformations has been performed by 
Grandati [I81 , but the full characterization of SL-OPS obtainable via multi-step Darboux transformations 
of mixed type remains an open problem. 

Another consequence of the conjecture is that all exceptional polynomials could be written as Wron- 
skian determinants involving essentially classical orthogonal polynomials (more specifically, involving one 
classical polynomial and many quasi-rational factorizing functions). 
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